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Abstract 

^ The dynamics of an initially localized wavepacket is studied for the generalized nonlinear 

^ Schrodinger Equation with a random potential, where the nonlinearity term is IV'I^V' ^md p is 

arbitrary. Mainly short times for which the numerical calculations can be performed accurately are 
^ considered. Long time calculations are presented as well. In particular the subdiffusive behavior 

^ where the average second moment of the wavepacket is of the form (7712) ~ is computed. Con- 

m 

^ trary to former heuristic arguments, no evidence for any critical behavior as function of p is found. 

""q,! The properties of a (t) are explored. 
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We consider the discrete nonhnear Schrodinger Equation (NLSE) with a random potential 
in one dimension: 

i-g^ = -i'n+l-'ipn-l+£ni'n + Pllpnflpn (l) 

where En are i.i.d. random variables uniformly distributed in the interval [— ^] , p is the 
degree of nonlinearity and (3 is its strength. For /? = this equation reduces to the Anderson 
model where all the states are exponentially locahzed [7]. Consequently, for /3 = if one 
starts with a localized wavepacket it will not spread indefinitely. In the absence of a random 
potential spreading takes place for all p |i4j. In fact the continuous version of ([T]) for p = 2 
and without the disorder is integrable [H]. The case of p = 2 is of experimental relevance in 
classical optics [H] and in the field of Bose-Einstein condensates, where the NLSE is known 
as the Gross-Pitaevskn equation [2l[l0]. This equation was studied extensively in the recent 
years, mainly for p = 2. In particular, the growth of the second moment was explored and 
it was found (numerically) to grow sub diffusively [U M, IS], namely, for a particle initially 
at n = 0, 

(m2 (t)) = Dr (2) 

where m2 = Xln^^ li^nf and a was found to be a ~ 0.33 (for p = 2). The average (...) 
denotes an average over the realizations of the random potential. 

The analytical and intuitive understanding of ([T| is quite poor. The simplest intuitive 
argument is that if a wavepacket spreads for long enough time, the amplitude of each state 
becomes negligible (since the norm, {ipuf = 1, is conserved) and as a result the nonlinear 
term weakens and becomes irrelevant, consequently, localization takes place. The difficulty 
with this argument (in addition with the fact that it disagrees with numerical results [H 
m [13]) is that although absolute value of the nonlinear term becomes smaller it should be 
compared to an energy scale that may decrease as well. For p = 2 such an argument was 
developed by Pikovsky and Shepelyansky [9] (It is very similar to an argument that was 
found to work remarkably well for another system |12j). We generalize this argument to 
an arbitrary value of p. Assuming that after some time the packet is spread over An 
states while its norm is preserved, than typically lipnf ~ ^ and therefore the nonlinear 
term produces an energy shift of the order 6E = (3 l^/'nT that is of the order of 5E ^ /\n'~^ . 
Comparing this term with the typical distance between the energies of the hnear problem, 
/\E gives ^ ^ (3/\n~^ . Based on this argument, Pikovsky and Shepelyansky that 
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were interested in the case p = 2, where ^ ^ (3 concluded that there is a critical value 
denoted by /3c such that for (3 < (3c the nonlinear term is negligible compared to the level 
spacings of the hnear problem and therefore Anderson localization holds. For (3 > I3c the 
levels of the linear problem are mixed and presumably Anderson localization breaks down 
and spreading takes place. From this argument it turns out that p = 2 is a critical degree 
of nonlinearity and for p > 2, ^ ^ as An grows and localization holds. Existence of a 
critical value of p was not considered in [H] since only the case p = 2 was studied. Also for the 
nonlinear Schrodinger equation without disorder, p = 2 has a critical meaning P,in]- In the 
present paper no evidence for the criticallity at p = 2 was found. This leads one to question 
the validity of the arguments implying the criticallity of p = 2 for spreading. Recently, 
Flach, Krimer and Skokos presented arguments that a = and there is no critical value 
of /? or p [H [5|, [13]. Their arguments are supported by some numerical calculations. It is 
unclear how a should behave when the limit p — > is taken since in this limit localization 
takes place and one expects a = 0. This is another motivation for the present work. Some 
arguments presented in [U [9l [I3] involve assumptions on chaoticity of various modes. The 
present work does not test these assumptions. 

There are conjectures based on perturbation theory [3| and rigorous results [15] claiming 
that asymptotically the second moment of the wave packet cannot grow faster than logarith- 
mically as a function of time. Nevertheless, numerical data predicts a power law growth of 
the second moment. If we trust the conjectures (their violation will be very surprising and 
of great interest) it is reasonable that the available numerical data is either not asymptotic 
(the time scale of this problem is unknown and therefore also the time when the system 
enters the asymptotic regime) or not reliable due to computational errors. Considering this, 
we concentrate on the short time behavior of a wave packet. Our results for the long time 
behavior are also presented for completeness. 

In order to follow the dynamics of a wave packet, we use the SABA algorithm, which 
belongs to the family of split step algorithms and evaluates the wave packet in small steps, 
changing from coordinate space to momentum space. We apply the disorder and nonlinear 
interaction in the coordinates space, transform the wave to momentum space and apply there 
the kinetic energy term, transform it back to the coordinate space and so on. Nearly all 
numerical calculations for this problem use such methods. Additional details on the SABA 
algorithm, can be found in reference [I3]. Like any numerical algorithm, the SABA algorithm 



3 



accumulates errors during the calculation which grow with the time of the integration. We 
use two criteria to determine whether our results are reliable or not: (tl) time reversal and 
(t2) comparison with data which is obtained using smaller time steps. Time reversal means 



end of this process (if there are no errors) we should get the initial wave packet. To measure 
the accumulated errors, we define 5tr = Ylm li'miuai — 'ipreversedl and demand 6tr < 0.1. The 
comparison with smaller time step is done as follows: we calculate the second moment m2 (t) 
for representative reahzations and then recalculate it using smaller time step (half of the 
original one). We define 



and demand < 0.01. Nearly all pubhshed numerical calculations used a more relaxed 
test: (t3) where in (t2) m2 is replaced by the average over realizations. 

We calculate m2 for various values of P and p and average over 5,000 realizations until 
time 1000. We verified that (tl) and (t2) are satisfied for representative reahzations. We 
use time steps of 0.1, 0.02, 0.01 and 0.00025 for /3 = 0.25, 0.5, 0.75 and 1, respectively, that 
were chosen to satisfy (tl) and (t2). In addition, data is presented for /3 = 2 and 4 using 
time steps of 0.1 where (tl) and (t2) are not satisfied. The results are shown in Fig. lb 
where a is obtained from fits similar to the one presented in Fig. la. Only the data in the 
interval 500 < t < 1000, that does not involve the initial spread was used in the fit of a. 
If we choose the time interval to be 300 < t < 1000 or 800 < t < 1000, our results do not 
change in a significant way. As we could expect, a {p ^ 0) ^0 and when p is large, a is 
very small. The maximal a is obtained for and nothing special happens for p = 2. We 
do not see any discontinuity for p = 0. All the lines in Fig. lb have similar shape and after 
forming linear transformations a = Cia + C2 where Ci and C2 are independent of a and p and 
depend only on /3, all the lines approximately coincide as shown in Fig. 2, leading us to the 
conclusion that there might be some scaling property. 

In our short time runs (t < 1000) the wavepacket didn't spread over many sites. In the 
case of maximal spreading (/3 = 4, p = ^) the second moment reached to a maximal value of 
150 and for the parameters P = 1, p = 2 the second moment was smaller then 55, while the 
localization length is about 6. When we follow the dynamics for longer times (for which (tl) 
and (t2) are not satisfied but (t3) is satisfied) the results support our previous conclusion 



integrating ([T]) from time to some later time and then integrating back to time 0. At the 
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Figure 1: (a) (7712) for /3 = 1 and p = 2 as a function of time. The blue solid curve is the second 
moment calculated numerically and the green dashed line is the fit which we use in order to find 
a. (b) a (p) for different values of /3. From top to the bottom: /? = 4, 2, 1, 0.75, 0.5, 0.25 (yellow 
stars, purple triangles, turquoise asterisks, red circles, green squares and blue diamonds). Only the 
data for /3 = 0.25, 0.5, 0.75 and 1 where points are connected by lines satisfy (tl) and (t2). For 
all realizations, W = 4 and maximal localization length is of 6 lattice sites. At the initial time the 
wavepacket populates one site (n = 0). 

that nothing critical happens for p = 2. We see that the wave packet spreads for all powers 
of nonlinearity p in a similar way, as shown in Fig. 3 for p = 0, 1.5, 2, 2.5, 4 and 8. Similar 
results are found in detailed studies of Mulansky [8]. 

In conclusion, we have found that for short times, there is no evidence of any critical 
phenomena neither for p = nor any p = 2. This conclusion is supported by long time 
calculations. In addition, we found that a has a maximum for p = | and there is evidence 
for scaling (Fig. 2). Understanding the physics of the a (p) plots, explaining why is the 
maximal spreading obtained for p = |, explaining of the origin of the scaling and finding 
the asymptotic behavior of a (p) are left for future research. 
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Figure 2: Fig. lb after rescaling. 
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